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A YIELD MANAGEMENT METHOD AND SYSTEM 

Field of Invention 

The present invention relates to a yield management 
5 method and system. 

Background of the Invention 

Yield Management Systems (YMSs) are commonly used as 
computer-based supports to facilitate decisions in 
several applications, such as in the transportation 
business and, in particular, in the airline industry. The 
goal of a yield management system is that of allocating 
an offered capacity to different types of requests 
(categories) , in order to optimise a yield parameter. For 
example, the yield management system is used to maximise 
the revenue that can be obtained from selling a given 
transportable weight and volume of a cargo flight; the 
objective is to fly an aircraft as profitable as possible 
without allowing low- fare freights to displace high- fare 
freights . 

All yield management systems use a forecast method 
to estimate the requests of the different categories; a 
decisional system is then applied to the estimated 
2 5 requests for determining a suggestion that should be 
given to maximise the yield parameter. Typically, 
stochastic or probabilistic methods are employed, in 
order to take into account an uncertainty of the 
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estimated requests (and not only their mean value) . 

A solution known in the art consists of estimating a 
simple probabilistic distribution of the variables 
5 involved in the yield management system. A stochastic 

linear programming model is then applied using the simple 
probabilistic distribution to determine an allocation 
model that maximizes a target function defining the yield 
parameter. An example of the aforementioned method is 
10 described in Williamson, E. L. "Airline Network Seat 
Inventory Control: Methodologies and Revenue Impact", 
Massachusetts Institute of Technology Phd Thesis - 
Cambridge (Mass) June 1992. 

15 A drawback of the solutions known in the art is that 

they become computationally intractable when a complexity 
of the problem defined by the yield management system 
increases. This is particular critical if the offered 
capacity is defined by more than one variable, such as 

2 0 the weight and the volume of a cargo flight. In these 

cases many simplifications, shortcuts and assumptions are 
introduced, in order to solve the problem in a reasonable 
amount of time. For example, the number of categories is 
limited or a manageable form of the probabilistic 

2 5 distributions is assumed (for example a Gauss, Poisson or 
uniform distribution) . Moreover, the known solutions 
discard some capacity variables (for example using only 
the weight and just checking a posteriori if volume 
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constraints are satisfied) or discard the stochastic 
nature of the yield parameter (for example considering 
only its mean value) . 

5 However, all the simplifications described above 

strongly reduce the accuracy of the yield management 
system. As a consequence, the solutions known in the art 
may lead to a relevant economic lost. 

10 It is an object of the present invention to overcome 

the above-mentioned drawbacks. In order to achieve this 
object, a method as set out in the first claim is 
proposed. 

15 Disclosure of the Invention 

Briefly, the present invention provides a yield 
management method for optimising a yield parameter 
resulting from assigning a capacity offered by a future 
instance of a service to each one of a plurality of 

20 different categories of requests competing for the 

capacity, the capacity being defined by at least one 
capacity variable, the method including the steps of: 
storing a set of historical profiles for each one of a 
plurality of previous instances of the service, the set 

25 including a historical profile of an historical value of 
each capacity variable reserved by each category, 
assigning a probability to each previous instance of the 
service, estimating a potential profile of a potential 
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value of the capacity variable from each historical 
profile according to a corresponding current value of the 
capacity variable reserved by the category for the future 
instance of the service and according to a corresponding 
5 unconstrained demand of the capacity variable for the 
category in the previous instance of the service, 
defining a historical scenario for each previous instance 
of the service, the historical scenario including a final 
potential capacity variable from each corresponding 
10 potential profile, and determining an authorisation to 

allocate the offered capacity for each capacity variable 
of each category in the future instance of the service by 
applying a stochastic model to the historical scenarios 
according to the corresponding probabilities. 

15 

Moreover, the present invention also provides a 
computer program application for performing the method, a 
program product storing the program application, and a 
corresponding yield management system. 

20 

Brief Description of the Drawings 

Further features and the advantages of the solution 
according to the present invention will be made clear by 
the following description of a preferred embodiment 
25 thereof, given purely by way of a non-restrictive 

indication, with reference to the attached figures, in 
which: 

Fig.l is basic block diagram of a data processing 
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system in which the yield management method of the 
invention can be used; 

Fig. 2 shows a partial content of a working memory 
included in the system; 

Figg.3a-3d are a flow chart describing the logic of 
the yield management method implemented in the system; 

Figg.4a-4b depict different diagrams of some 
variables involved in the method. 

Description of the Preferred Embodiments 

With reference in particular to Fig.l, there is 
shown a data processing system 100, which is used as a 
computerised support for controlling cargo flights of an 
airline. The system 100 is formed by a reservation 
subsystem 105 and an authorisation subsystem 110; the 
subsystems 105 and 110 are connected to each other 
through a communication channel 115. 

The reservation subsystem 105 consists of a 
mainframe having a Central Processing Unit (CPU) 120. A 
working memory (RAM) 125 and a bulk memory 130 (generally 
consisting of one or more hard-disks) are associated with 
the CPU 120. Moreover, a telecommunication control unit 

(TCU) 135 is used to connect multiple remote terminals 
140 and the communication channel 115 to the CPU 120. The 
reservation subsystem 10 5 hosts a transactional 
processing environment, which allows the remote terminals 

140 to respond immediately to each request entered by an 
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The authorisation subsystem 110 consists of a 
workstation having a RISC architecture. The authorisation 
subsystem 110 includes a Central Processing Unit (CPU) 
150 and a working memory (RAM) 155 that is used directly 
by the CPU 150. The CPU 150 is associated with a bulk 
memory consisting of a hard-disk 160 and of a driver unit 
165 for reading CD-ROMs 170. Moreover, the CPU 150 is 
connected to an input unit (IN) 175, which consists for 
example of a keyboard and a mouse, and to an output unit 
(OUT) 180, which consists for example of a monitor and a 
printer. A network interface card (NIC) 185 is used to 
connect the communication channel 115 to the CPU 150. 

The above-described architecture implements a yield 
management system. Yield management systems are a well- 
known application field of several Operation Research 
(OR) techniques. Basically, the goal of a yield 
management system is that of optimising one or more 
parameters that can be obtained from a given capacity of 
a generic service. A demand consisting of multiple types 
of requests (categories or classes) competes for the 
offered capacity; each category is defined by a subset of 
the requests exhibiting similar characteristics (such as 
a physical property, a fare, a reservation/cancellation 
behaviour) . 
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The yield management system selects the requests by 
accepting only the more convenient ones; this is done by 
assigning a booking limit (authorisation) to each 
category. The offered capacity may be allocated to the 
5 different categories using several techniques. In an out 
of nesting policy, the offered capacity is simply split 
among the categories according to their authorisations. 
This is a very conservative technique, which ensures a 
pre-set capacity to each category; however, it may lead 

10 to a yield loss when the requests for a category are 

lower than the respective authorisation. Conversely, in a 
nesting policy the same capacity may be reserved by more 
than one category; the categories are ranked according to 
their value, in order to assign a greater capacity to the 

15 most valuable ones. Particularly, in a total nesting 
policy each category may reserve all the capacity 
assigned to the category itself and to the categories 
with a lower value. This is a very aggressive technique, 
which favours the most valuable categories to the 

2 0 detriment of the other ones; however, it provides good 

results only when the requests for low- fare freights are 
received before the ones for high- fare freights. 
Otherwise, in a partial nesting policy each category may 
reserve the capacity assigned to the categories with a 

2 5 lower value only when the capacity assigned to the 
category itself has been completely allocated. 



In the example at issue, the yield management system 
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is used to maximise a revenue A that can be obtained 
from a given capacity that is offered by each cargo 
flight; the offered capacity is defined by two different 
variables consisting of a total weight P tot and of a total 
5 volume v to t of the freights that may be transported by the 
flight. The freights are classified in five categories Ci 
(with i=1...5), for example defined as normal light, 
special light, normal heavy, special heavy, and extra; 
this number is a good compromise between the opposed 
10 requirements of having uniform characteristics within 
each category (leading to a high number of categories) 
and a constant behaviour of each category (leading to a 
low number of categories) . 

15 The yield management system employs a partial 

nesting policy, with two sets of independent 
authorisations. A set of nesting authorisations Ap ± is 
used for the weight, with the categories that are ranked 
in a respective descending nesting order ORDp (from 1 to 

20 5); a set of nesting authorisations Avi is used for the 
volume, with the categories that are ranked in another 
independent nesting order ORDv. The partial nesting 
policy is preferred in this application since the period 
on which each flight is open for reservation is short, so 

25 that the order of the requests is very random. Moreover, 
it is generally difficult to rank the freights according 
to their value in an accurate manner, because both the 
weight and the volume are to be taken into account. In 



GB920010034US1 9 

addition, the weight and the volume are not always 
constraints limiting the offered capacity at the same 
time . 

5 As described in the following, the weight nesting 

authorisations Ap i# the volume nesting authorisations 
Avi, the weight nesting order ORDp and the volume nesting 
order ORDv for a future instance of each flight are 
determined by the authorisation subsystem 110. The 

10 authorisation subsystem 110 further supplies a minimum 
selling price, which consists of a couple of values 
indicating a weight unitary revenue MAp and a volume 
unitary revenue MAv for each segment of a network 
supported by the airline, under which a request should be 

15 refused because not profitable. The nesting 

authorisations Api,Avi, the nesting orders ORDp, ORDv and 
the minimum selling price MAp,MAv are transmitted to the 
reservation subsystem 105 and stored on its hard disk 
130. The hard disk 130 further stores booking information 

2 0 BK about every request accepted by the reservation 

subsystem 105 for the future instance of each flight, 
such as a date of the accepted request, a weight and a 
volume of the freights, a name of the client, an applied 
fare (if already available) . The sum of the weight and 

25 volume of the accepted requests gives a weight Rpi and a 
volume Rvi, respectively, already reserved by each i-th 
category. The reservation subsystem 105 calculates a 
residual weight Dpi and a residual volume Dvi that are 
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still available for the i-th category according to the 
following formulas: 



The operator at a remote terminal 140 enters a 
reservation request for a particular flight. The 
reservation request is accepted if both the weight and 
the volume of the freights do not exceed the residual 
weight Dpi and the residual volume Dvi, respectively, of 
the corresponding category; otherwise, the reservation 
request is denied. At the same time, the reservation 
subsystem 105 adds together the minimum selling price of 
all the segments from an origin to a destination of the 
request, and multiplies the result by the weight and the 
volume of the freights; when the revenue of the request 
(if available) is lower than the value so obtained, the 
request is always refused. When the reservation request 
is accepted, it is recorded on the hard disk 130 and the 
respective residual capacities Dp if Dvi are accordingly 
updated. The booking information BK, the reserved 
capacities Rpi,Rvi and the residual capacities Dpi,Dvi for 
the future instance of each flight are sent to the 
authorisation subsystem 110 either periodically (for 
example at the end of each day) or on request. 
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Similar considerations apply if the data processing 
system has a different architecture, for example 
consisting of a single computer or with the reservation 
subsystem and the authorisation subsystem working in the 
5 INTERNET, if the reservation subsystem and the 

authorisation subsystem have a different structure, for 
example if they consist of Personal Computers (PCs) or 
they are implemented with a local area network, if the 
freights are classified in a different number of 
10 categories, if other categories are envisaged, if a 

different nesting policy is employed (for example a mixed 
partial one) , and so on. 

Considering now Fig. 2, there is shown a partial 
15 content of the working memory 155 of the authorisation 
subsystem; the information (programs and data) is 
typically stored on the hard-disk and loaded (at least 
partially) into the working memory 155 when the programs 
are running, together with an operating system and other 
2 0 application programs (not shown in the figure) . The 

programs are initially installed onto the hard disk of 
the authorisation subsystem from CD-ROM. 

An input /output interface (I/O) 205 controls 
25 interaction of a space controller with the authorisation 
subsystem. A network interface (NET) 210 processes a set 
of protocol layers, which allows the authorisation 
subsystem to exchange information with the reservation 
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subsystem (through the communication channel) . 

The working memory 155 stores a historical 
repository 215 for each flight. The historical repository 
5 215 logs information about a series of previous instances 
of the flight; preferably, the last 104 previous 
instances of the flight (one per week in two years) are 
logged. Particularly, the historical repository 215 
includes a field INFOj (with j=1...104) that identifies the 

10 previous instance of the flight (for example its number 
and departure time) , and a field SPHj with information 
about each shipment (for example the weight and the 
volume of the freights that have been transported, and a 
fare actually applied) . The historical repository 215 

15 further stores a set of historical profiles, which 

consist of a temporal profile of the weight and of the 
volume that have been reserved by each category during an 
opening period of the previous instance of the flight. 
Each historical weight profile and each historical volume 

20 profile is defined by a series of daily snapshots (for 
example taken at the end of the day) ; if the previous 
instance of the flight was open 2 5 days before its 
departure, the k-th snapshot (with k=-2 5...0) for the i-th 
category in the j-th historical weight profile consists 

2 5 of a reserved weight Rpi jk and a residual weight Dp ijk , 

whereas the same snapshot in the j-th historical volume 
profile consists of a reserved volume Rvij k and a residual 
volume Dv ijk ; the first snapshot gives the reserved 
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capacity RPij-25/ RVj.j-25 and the residual capacity Dpij_ 
25 ,Dvij_25 at the opening of the previous instance of the 
flight, whereas a last snapshot gives the reserved 
capacity RpijO/ Rv ij0 and the residual capacity Dpij 0 ,Dvij 0 at 
5 its departure. 

The working memory 155 further stores a current 
repository 22 0 for the future instance of each flight. 
Particularly, the current repository 220 includes a field 

10 INFO that describes the future instance of the flight 
(for example the offered capacity P t ot,V tot and its 
scheduled departure time) , a field with the corresponding 
booking information BK, and a field with a current time t 
at which the information is taken. The current repository 

15 220 further stores a set of current profiles, which 

consist of a temporal profile of the weight and of the 
volume that have been currently reserved by each category 
for the future instance of the flight. Each current 
weight profile and each current volume profile is defined 

20 by a series of daily snapshots. The k-snapshot (with k 
ranging from -2 5 to the current time t) for the i-th 
category in the current weight profile consists of a 
reserved weight Rpi k and a residual weight Dp ik , whereas 
the same snapshot in the current volume profile consists 

25 of a reserved volume Rv ik and a residual volume Dv ik ; the 
last snapshot gives the reserved capacity Rpi t ,Rv it and 
the residual capacity Dp it ,Dvi t at the current time t. 
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The booking information BK, the current reserved 
capacities Rpi t ,Rv it and the current residual capacities 
DPit,Dv it are received from the reservation subsystem 
through the network interface 210. The set of profiles 
5 stored in the current repository 220 for each flight are 
copied to the corresponding set of profiles in the 
historical repository 215 after departure of the flight; 
the shipment field SHPj is continually updated by the 
input/output interface 2 05 as soon as the corresponding 
10 information becomes available. 

The historical repository 215 and the current 
repository 220 are accessed by a weighting module 225, 
which assigns a weight Wj to each previous instance of 
15 the flight. The weight Wj depends on a relevance of the 
previous instance of the flight with respect to its 
future instance. 

The weights Wj are input to an estimation module 
20 23 0, which further accesses the historical repository 215 
and the current repository 220. The estimation module 230 
determines a potential (weight and volume) profile from 
each corresponding historical profile. The potential 
profile is obtained adjusting the historical profile 
25 according to the respective current profile and an 

estimated unconstrained (or untruncated) demand of the 
respective capacity variable for the category in the old 
instance of the flight. The unconstrained demand takes 
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into account the periods of time in which the category 
has been closed to sales (because the reserved weight 
and/or the reserved volume for the category have reached 
their booking limits defined by the respective nesting 
authorisations) , and estimates a demand not limited by 
the offered capacity. The estimation of the unconstrained 
demand is carried out independently for the weight and 
for the volume. Each potential profile consists of 25 
snapshots; the k-th snapshot for the i-th category in the 
j-th potential weight profile consists of a potential 
weight request Rp ijk , whereas the same snapshot in the j - 
th potential volume profile consists of a potential 
volume request Rv ijk . The last snapshot of the potential 
weight profile and of the potential volume profile define 
a final potential weight request Rp iJ0 and a final 
potential volume request Rv, j0 , respectively, for each i- 
th category in the j-th previous instance of the flight. 

The estimation module 23 0 further determines a 
potential revenue A~ for each i-th category transported 
by the j-th previous instance of the flight. The 
potential revenue is estimated according to the revenue 
of the category in the previous instance of the flight; 
this value is further adjusted to take into account the 
revenue of the corresponding capacity currently reserved 
for the future instance of the flight. 
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The estimation module 23 0 outputs a series of 
historical scenarios for each flight. The j-th scenario 
consists of the corresponding final potential weight 
request Rp lJ0 , final potential volume request Rv.jo , and 
potential revenue A~ for each i-th category. The space 
controller defines a user scenario through the 
input/output interface 205. The user scenario consists of 
a potential weight request Rp i , a potential volume 
request Rvt , and a potential revenue A ; . , which are 
independently forecast by the space controller for each 
i-th category. A weight w is assigned to the user 
scenario by the space controller through the input/output 
interface 205. 

The weight Wj of each historical scenario 
{Rp ij0 ,Rvuo and the weight w of the user scenario 

{Bp l , Rvi , A^}, which will be denoted as a whole with w s in 
the following, are input to a normalisation module 235. 
The normalising module 235 determines a probability w s 
of each (historical and user) scenario. 

The probabilities w s are input to a nesting order 
module 240, which also receives the potential revenues 
XT of each historical scenario from the estimation 
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module 23 0 and the potential revenues A,, of the user 
scenario from the input/output interface 205. An input 
weight nesting order IN_ORDp and an input volume nesting 
order IN_ORDv are provided by the space controller 
thorough the input /output interface 2 05, and then 
transmitted to the nesting order module 240. The nesting 
order module 240 outputs the nesting orders ORDp,ORDv for 
the future instance of the flight. 

The nesting orders ORDp,ORDv are provided to an 
optimisation module 245; the optimisation module 245 
further receives the historical scenarios { Rp ij0 , Rv,jo , A y . } 
from the estimation module 23 0 and the user scenario 
{Rp i , R~Vi , X7} from the input/output interface 205, 
together with the respective probabilities w s from the 
normalisation module 235. The space controller supplies 
an aggressiveness parameter AG, indicative of an attitude 
to the risk, to the optimisation module 245 through the 
input/output interface 205. The optimisation module 245 
employs stochastic programming techniques to determine 
the (optimal) nesting authorisations Api,Avi that 
maximise a potential revenue A of the future instance of 
the flight. 

Stochastic programming techniques are commonly used 
as a support in decision-making processes under 
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uncertainty, such as in financial trading. An example of 
stochastic programming technique is described in "IBM 
Stochastic Solution", Tito A. Ciriani, Stefano Gliozzi, 
AIROnews Bulletin of AIRO, VI, n.l - Spring 2000, pp. 6- 
10. In particular, the module 245 optimises each scenario 
using a deterministic linear programming technique; the 
results are then combined weighting each scenario 
according to the respective probability. 

Similar considerations apply if a whole application 
program (defined by the different modules described above) 
and the corresponding data are structured in a different 
manner, if other functions are provided, if the historical 
and current repository are replaced by equivalent memory 
structures, and so on. 

With reference now to Figg.3a-3d, the application 
program when running on the authorisation subsystem 
performs a method 300. The method starts at block 302 and 
then passes to block 304, wherein a menu with a series of 
possible choices is displayed on the monitor of the 
authorisation subsystem. 

The method carries out the operations corresponding 
to the selected choice. Particularly, if the space 
controller has selected a management function, the blocks 
3 05-307 are executed, whereas if the space controller has 
selected an optimisation function, the blocks 308-374 are 
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executed; conversely, if the space controller has 
selected an exit option the method ends at the final 
block 376. 



5 Considering now block 305 (management function) , the 

information about the future instance of the flight is 
received from the reservation subsystem, and it is stored 
into the current repository. The method passes to block 
3 06, wherein the set of profiles stored in the current 

10 repository are copied to the historical repository after 
departure of the flight. Continuing to block 307, the 
information about a selected shipment for the previous 
instance of the flight is updated as soon as it is 
available. The method then returns to block 3 04 waiting 

15 for a new command. 



With reference to block 3 08 (optimisation function) , 
a flight is selected either automatically (in a batch 
processing of all the flights of the airline) or by the 
2 0 space controller (in an interactive processing) . 



The method then assigns the corresponding weight Wj 
to the previous instances of the flight. Particularly, a 
season coefficient wsj for each j-th previous instance of 
25 the flight is calculated at block 308. The coefficient 
wsj depends on a season of the previous instance of the 
flight, that is on a difference between a period of the 
previous instance of the flight and its future instance; 
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the more the previous instance of the flight has happened 
in a similar period of the year as its future instance, 
the more the season coefficient wsj is bigger. The season 
coefficient wsj is obtained from a matrix having 52 rows 
5 and 52 columns, which is preferably symmetric; each 

element [u,v] of the matrix (with u,v=1...52) stores the 
season coefficient wsj of a previous instance of the 
flight that departed in the u-th week of the year, for a 
future instance of the flight that will depart in the v- 
10 th week of the year; the season coefficients wsj have a 

value ranging from 0 to 1, with all the diagonal elements 
of the matrix (u=v) that are equal to 1. 

Descending into block 312, a time distance 
15 coefficient wgj for each j-th previous instance of the 
flight is calculated. The coefficient wgj depends on a 
space between the current time and the departure of the 
previous instance of the flight; the more recent the 
previous instance of the flight is, the bigger the time 
20 distance coefficient wgj is. The time distance 

coefficient wgj is defined as a hyperbolic function of a 
number of days g that have passed from the departure of 
the previous instance of the flight. We set wgj=l for g=0 
and wgj=ai for g—>°° (with 0<ai<l, and preferably ai=0) ; 
2 5 moreover, we set an exponent a 2 of the exponential law 
(with a 2 <l, and preferably a 2 =0.2), and a scale 
coefficient a3 for reducing a concavity of the hyperbola 
(with a 3 >30, and preferably a 3 =30) . The time distance 
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coefficient wg 3 - is then given by: 



1 



5 




This results in a time distance coefficient wg-j=0.53 
for g=708 (2 years) , so as to give more relevance to the 
most recent previous instances of the flight, while 
10 maintaining some relevance for the oldest previous 
instances of the flight at the same time. 

The weight Wj for each j-th previous instance of the 
flight is then calculated at block 314, by combining the 
15 corresponding season coefficient wsj and time distance 

coefficient wg-j . Particularly, the value Wj is calculated 
as a sum of the coefficients wSj and wgj , which are 
weighted according to a parameter h with a value between 
0 and 1, and preferably equal to 0.65: 

20 

Wj = h ■ WSj +(l-h)- wgj 

Considering now block 316, an opening coefficient 
aijk for each k-th snapshot of the i-th category in the j- 
2 5 th previous instance of the flight is estimated; the 

opening coefficient a i; j k indicates a percentage of time 
during which the category has been open to sales during a 
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period associated with the snapshot, that is from a 
preceding snapshot (k-1) to a current snapshot (k) . 

As shown in Fig. 4a, the category is open only when 
5 the residual weight and the residual volume are both 
positive. We denoted with s + (f) a function that 
determines the sign of its argument f, that is: 



The starting conditions that are used (in the given 
order) to determine whether the category has been always 
open or always closed are: 



Otherwise, we calculate the sign of the derivative 
of the functions defined by the residual weight and the 
residual volume: 
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1 if ^(Dp gk _ l ) = s\Dp. k ) = s + (DVy*J = * + (Dv ijk ) = 1 

0 if s + {Dp ljk _ x ) = s + (D PiJk ) = 0ors + {Dv ijk _, ) = s + (Dv ljk ) = 0 



20 



sp\j k = s+ ( D P ijk )~ s+ ( D P tjk -i ) 
sv' ijk =s + (Dv ijk )-s + (Dv lJk _ 1 ) 



When sp' ijk ,sv' ijk *0 , the intersection of these functions 
with a time axis is given by: 
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Fljk Dp i]k _,-Dp ljk 



We now define the following functions: 



\-tp ijk if sp' ijk =+\ 
tv iJk if sv' iJk = -1 



If we consider all the possible cases, we can 
demonstrate that: 

15 

tPijk if Wok * 0 and sv' ijk = 0 
K if sv ijk * 0 and sp iJk = 0 

min-^V, tv ijk }ff sp'et ' sv \jk > 0 
maxjo, (tp ijk + tv iJk ) - l\otherwise 



The method continues to block 318, wherein a weight 
emphasis value Ep ik and a volume emphasis value Ev ik to be 
applied to the k-th snapshot for each i-th category 
25 (irrespective of the old instance of the flight) are 

calculated; the emphasis values Epi k ,EVi k indicate what 



GB920010034US1 



24 



could have happened if the category had been always open 
to sales in the period associated with the k-th snapshot. 
The emphasis is calculated as a function of observations 
consisting of the demand of the respective category in 
5 the corresponding historical profiles where the category 
was open in the same period. 

The method takes into account all the corresponding 
(weight or volume) historical profiles for the same 

10 category, which were always open in the period (aijk=l) . 
Further historical profiles (with a decreasing value of 
the opening coefficient aijk) are also taken into account, 
if necessary, in order to ensure a sufficient number of 
observations (for example at least equal to 30% of the 

15 number of the previous instances of the flight) . If we 
denote with P' l]k = P ijk - P ijk _ x and with V ljk =V ijk -V ljk _ { a gradient 
of the reserved weight and of the reserved volume, 
respectively, the emphasis values Epi k ,Evi k are defined as 
a weighted mean of these values : 

20 

Ep*=^ 

j 

£v*=-A^ 

2 5 2w*V 

j 
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Passing to block 32 0, the method estimates a 
potential weight gradient Upij k and a potential volume 
gradient Uvijk for each k-th snapshot of the i-th category 
in the j-th historical weight profile and historical 
5 volume profile, respectively. The potential gradients 

Upijk/Uvijk are obtained as a linear interpolation between 
the value Py k ,V yk for aijk=l and the value 
max[Py k ,Ep ljk \max^ jk ,Ev. k ] for a ijk =0, that is: 

Up iJk = max[i^ , E PiJk + a iJk + (P yk - Ep ljk )] 
Uv iJk = max[7 # , Ev iJk + a ljk + (V' lJk - Ev ljk )] 



15 The method then builds each potential profile. As 

shown in Fig. 4b, the potential profile is set to the 
respective current profile up to a time t' corresponding 
to the current time t. Considering now block 322, the 
next snapshot for the i-th category in the j-th potential 

20 weight profile is calculated integrating the potential 
weight gradient Upijk, whereas the same snapshot for the 
i-th category in the j-th volume profile is calculated 
integrating the potential volume gradient Uvij k ; the 
integration is carried out with a simple sum, that is 

2 5 Rp ijk = Rpy^+Up^ and Rv y k = RVi/k-i + Uv ijlc , respectively. The 

potential gradients relating to the period including the 
time t ' are reduced to take into account the time that 
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has passed from the preceding snapshot; we denote with k 
an index of the snapshot immediately following the time 
t'; if g t and g- are the number of hours from the time t' 

and from the k-th snapshot to the departure of the 
5 previous instance of the flight, respectively, the 

potential gradients L^-,£/v.- are reduced according to the 
following formulas: 



The method then checks at block 324 whether either 
the resulting potential weight request P yk or the 
resulting potential volume request V ijk for the current 

15 snapshot is not strictly positive (>0). If so, the method 
descends into block 326, wherein the potential weight 
request P ljk and the potential volume request V ijk are both 
set to zero; the method then continues to block 328. 
Conversely, if both the resulting potential weight 

2 0 request P iJk and the resulting potential volume request 

V ijk are strictly positive, the method passes to block 328 
directly. In other words, the potential weight capacity 
P iJk and the potential volume capacity V ljk (with k=t,&,...0) 
for the current snapshot are given by: 
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Pgk-i + Up i]k if Pijk-i + Up ijk > 0 and V ijk -\ + Uv ljk > 0 
0 otherwise 



Vijk-i + Uv ijk if Pijk-i + Up ijk > 0 and V, Jk -i + Uv iJk > 0 
0 otherwise 



With reference now to block 32 8, the result is 
further corrected by reconciling the potential weight 

request P Ijk and the potential volume request V ijk of the 
10 current snapshot, in order to resolve any inconsistency 

that may be introduced between the two capacity variables 
and that may bring about incongruent values of a density 

p uk 

of the potential request. We denote with 0^=—^ the 

V i jk 

density of the reserved capacity for the i-th category in 
15 the k-th snapshot of the j-th previous instance of the 

flight, and with Syk =-^L the density of the potential 

Vijk 

capacity for the same snapshot. If the density Syk for 
the current snapshot is comprised between 5ijk and 8i }k -\ , 
not reconciliation is carried out and the method descends 
20 into block 33 0 (described in the following) . 



Otherwise, the potential demand is correct to an 
admissible reference density R8ij k at block 332. The 
reference density R8ij k is set to the value between 8ij k 
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and Sijk-i that is closer to 8yk . Particularly, the 
potential weight request Rp ljk and the potential volume 
request Rp ijk are reconciled according to the following 
formula : 

5 p ijk ^l(p. k+ v ijk .R8 ijk ) 

Considering now block 33 0, the method checks whether 
10 the last snapshot (associated with the departure of the 

old instance of the flight) has been reached. If not, the 
method returns to block 324 for processing a next 
snapshot . 

15 Conversely, the method descends into block 333, 

wherein a historical revenue Aij for the i-th category 
transported by the j-th previous instance of the flight 
is calculated by summing a value of the respective 
shipments . The historical revenue Aij may be calculated 

20 only when the value of all the shipments for the category 
are known; otherwise, if the value of one or more 
shipments is not available, the historical revenue Aij 
for the whole category is deemed unknown. In a similar 
manner, a current revenue Ai for the i-th category in the 

2 5 future instance of the flight is calculated by summing a 
value of the respective requests accepted at the current 
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time t. 

In the cargo flights, the freight fares are defined 
according to a so-called taxable weight P*; the taxable 
5 weight P* is a virtual weight of the freights that is 

applied to calculate their value for a given unit weight 
fare. We denote with A* a reference density (for example 
A*=167 kg/m 3 ) . If a density 8 of the freights is equal to 
or higher than A* (heavy freights) the taxable weight P* 
10 of the freights is equal to their real weight; otherwise, 
if the density 5 of the freights is lower than A* (light 
freights) , the taxable weight P* of the freights is 
obtained multiplying their volume by the reference 
density A*: 

15 

P*=max(P,V-A*) 

A (weight) unit revenue X is then calculated as X = -^r 

P 

We define a historical taxable weight for each i-th 
2 0 category in the j-th previous instance of the flight as 
Rp* 0 =m.ax(Rp lJ0 ,Rv ij0 ■ A*), and a current taxable weight for 
each i-th category as Rp* u = max(flp f , , Rv u ■ A* ) ; it is then 
possible to calculate a historical unit revenue Adj for 
each i-th category in the j-th previous instance of the 
2 5 flight and a current unit revenue X± for each i-th 

category in the future instance of the flight by applying 
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the above mentioned formula, that is X- = — |- and 

A t = — — , respectively . 
R Pu 



Proceeding to block 334, the historical unit 
5 revenues X±j , for each i-th category in the j-th previous 
instance of the flight wherein the respective historical 
revenue Aij is unknown, are set to the weighted mean of 
the respective values Xij available, that is: 




The method then estimates a potential unit revenue 
X y for each i-th category in the j-th previous instance 

15 of the flight, by adjusting the corresponding historical 
unit revenue A,ij according to the corresponding current 
unit revenue X±. Particularly, a departure coefficient wd 
for the future instance of the flight is calculated at 
block 33 6. The coefficient wd depends on a space between 

2 0 the current time t and the departure of the future 

instance of the flight. The more the current time is 
close to the departure, the bigger the coefficient wd is. 
We assume that the current unit revenue Xi does not 
affect the historical unit revenue X±j when the number of 



GB920010034US1 



31 



days to the departure of the future instance of the 
flight is greater than a pre-set value g m (for example 
g m =10) . If we denote with g d the number of days from the 
current time t to the departure of the future instance of 
5 the flight, the departure coefficient wd is given by: 



Passing to block 338, an increment coefficient wxij 
for each i-th category in the j-th previous instance of 
the flight is further calculated. The coefficient wxij 
depends on an increment of the potential capacity at the 
departure for the i-th category in the j-th previous 
instance of the flight with respect to the current 
reserved capacity. We define a potential taxable weight 
for each i-th category in the j-th previous instance of 

the flight as P* 0 = max(Rp IJ0 ,Rp ljO • A*); the increment 
coefficient wx^ is then given by the following formula: 



0, '4=; " 

pL 



0 otherwise 



ifp; 0 >o 



A corrective factor wcij for each i-th category in 
25 the j-th previous instance of the flight is then 

calculated at block 340, by combining the corresponding 
departure coefficient wd and increment coefficient wxij . 
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Particularly, the corrective factor wcj is calculated as 
a sum of the coefficients wd and wxij, which are weighted 
according to a parameter he having a value between 0 . 5 
and 1, and preferably equal to 0.6: 
5 wcij=hc-wxij+ (1-hc) -wd 

The potential unit revenue X y for the i-th category 
in the j-th old instance of the flight is then calculated 
at block 342 as a sum of the corresponding historical 
10 unit revenue A-ij and current unit revenue X±, which are 
weighted according to the corrective factor wcij , that 
is : 

\ = wc ;/ • Xy + (1 - wc iy ) • \ 

15 

The potential revenue A SJ for each i-th category in 
the j-th old instance of the flight is simply obtained as 
Ay = Xg • Rp* j0 . 

20 

Proceeding to block 344, the historical scenarios 
{Rp lj0 , Rvijo , A v } are determined for each old instance of 
the flight. The method then checks at block 348 whether 
the space controller wishes to input a user scenario. If 
25 not, the method descends into block 350 (described in the 
following) . Conversely, the space controller inputs the 
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user scenario {Rp, f Rvi , A, } at block 352 and inputs its 
weight w at block 354; the method then passes to block 
350. 

Considering now block 350, the probability w s of 
each {historical and user) scenario is calculated 
normalising the corresponding weight w s , that is: 



2>: 

10 The method then proceeds to block 355, wherein the 

input nesting orders IN_ORDp, IN_ORDv are provided by the 
space controller. Descending into block 3 56, the method 
checks whether for each i-th category some of the j-th 
historical scenarios have the final potential weight 

15 request Rp. J0 , the final potential volume request Rv lj0 and 
the potential revenue A tj all strictly positive. If not, 
the method passes to bock 3 58 (described in the 
following) . Conversely, the method enters block 3 60, 

wherein a potential weight revenue hp. = , and a 

2 0 potential volume revenue Av & ==^=r are calculated for 

each i-th category that meets the above mentioned 
conditions (taking into consideration all the historical 
scenarios and the user scenario) . A corresponding mean 
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weight revenue 7q) i and a corresponding mean volume 

revenue Av ; . are then calculated for these categories as 
a weighted mean of the potential weight revenues 7<p is and 
the potential volume revenues Xv is , respectively: 



The method then passes to block 3 58. 

Considering now block 3 58, the position of each 
category in the nesting orders ORDp,ORDv is then 
determined. For each capacity variable (weight and 
volume) , the categories that do not have a respective 
mean revenue maintain their position in the input nesting 
orders IN_ORDp , IN_ORDv . The other categories having a 

respective mean revenue Ap t , /Lv, are ranked in a 
corresponding descending order within the available 
positions. For example, let us suppose that the input 
weight nesting order for the five categories is 



4ft 
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IN_ORDp={C 2 , C 4 ,C 3 ,C 5 ,Ci}. The categories C 2 , C 3 and C 5 are 
ranked in the order {C 3 ,C5,C 2 } according to their 

respective mean weight revenues , Ap 5 , 7q> x ; conversely, 

the mean weight revenue Ap 4 (for the category C 4 ) and the 

mean weight revenue Ap { (for the category Ci) are not 
available. In this case, the resulting weight nesting 
order is ORDp={C 3 , C 4 , C 5 , C 2 , C x } . 

The method then proceeds to block 3 62, wherein the 
aggressiveness parameter AG is input by the space 
controller. The method builds a stochastic model that 
maximises a target function defined by the potential 
revenue A of the future instance of the flight; the 
potential revenue A is calculated using both a nesting 
policy (irrespective of an order of the requests) and an 
out of nesting policy; the nesting policy (best case) and 
the out of nesting policy (worst case) are combined 
according to the aggressiveness parameter. 

Particularly, a portion of the target function for 
the nesting policy is built at block 364. We define a 

potential nesting constrained weight TNp^ and a 

potential nesting constrained volume TNv i that would 
have been transported for each i-th category by the j-th 
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previous instance of the flight complying with the 
corresponding nesting authorisations Api and Avi, 
respectively: 

5 * , 

^TNv y <Av, 

<i 

The variables TNp gf TNv y are subjected to the 
• TNv„ 

constraint TNp ij =-=^- and to the simple bounds 

S ijO 

10 j^<p off ,7m:<v off . 

A potential nesting revenue TNAj of the j-th 
previous instance of the flight that would have been 
produced by the potential nesting constrained weight 

15 TNp y . and the potential nesting constrained volume TNv tJ 
is given by: 

The method then builds a further portion of the 
2 0 target function for the out of nesting policy, using two 
independent components for the weight and for the volume. 
Particularly, a set of out of nesting weight 
authorisations AOpi and a set of out of nesting volume 
authorisations AOvi are defined at block 366. If we 
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10 



denote with Pe the offered weight that is not used, and 
with Ve the offered volume that is not used, the out of 
nesting authorisations AOpi,AOvi are defined by the 
following relations: 

P = ^AO Pi +Pe 

5 

V = y ^AOv i +Ve 

/=i 

The out of nesting authorisations AOpi and AOvi are 
further constrained by the nesting authorisations Api and 
Avi, respectively: 



15 ^Ap t =^AOp t +Pe 

1=1 Si 

5 

Y, Av i = ^AOv i +Ve 

i=i <i 

We now define a potential out of nesting constrained 
2 0 weight TOp l} that would have been transported for each i- 
th category by the j-th previous instance of the flight 
in the worst case complying with the out of nesting 
authorisations AOpi,AOvi only for the weight, that is: 
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TOpy < AO Pi 
7Vp~:<AOv l -I~l 

Even in this case, the potential out of nesting 
constrained weight TOp tJ is subjected to the simple 
5 bounds T0p y <P off . A potential out of nesting weight 
revenue TOApj of the j-th previous instance of the 
flight that would have been produced by the potential out 
of nesting constrained weight Top tj is given by: 

10 TOApj^fdp^.^ 

!=1 

Proceeding to block 3 68, we now define a potential 
out of nesting constrained volume TOv tJ that would have 
15 been transported for each i-th category by the j-th 
previous instance of the flight in the worst case 
complying with the out of nesting authorisations 
AOpi,AOvi only for the volume, that is: 

TOv i} < AOv, 

As above, the potential out of nesting constrained 
volume TOv y is subjected to the simple bounds TOv.. < V lot . 



GB920010034US1 



39 



A potential out of nesting volume revenue TOAvj of the 
j-th previous instance of the flight that would have been 
produced by the potential out of nesting constrained 

volume Tov tJ is given by: 



The method then passes to block 370, wherein the 
target function to be optimised is defined combining the 
nesting portion and the weight and volume component of 
the out of nesting portion with the respective 
probabilities, according to the aggressiveness parameter 
AG: 



The nesting authorisations Api,Avi that maximise the 
target function defined above (with the variables meeting 
the respecting constrains and single bounds) are then 
obtained at block 372. 

Several by-products of this process may also be 
obtained. For example, the stochastic model calculates 
the minimum selling price MAp,MXv; particularly, the 
value MXp is equal to a marginal weight revenue and the 



TOAv J =f,TOv iJ -Zp r S i 




-) • 2(r0Ap, • w s ) + (1 -4^) • ^(TOAv s ■ w s ) 
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value MAv is equal to a marginal volume revenue in the 
neighbourhoods of the optimal solution, which are 
obtained as the partial derivative of the target function 
with respect to the offered weight and offered volume, 
5 respectively. Moreover, the stochastic model calculates a 
mean potential constrained weight CRp y . and a mean 
potential constrained volume CRv y as: 

CRp i} = AG-^(TNp > s ■w~ s ) + (l-AG) y £(TOp s -wj 
CRV. = AG • V (JWv, • V s ) + (1 - AG) • Y (ToV s ■ V s ) 

10 s 

The method verifies at block 374 whether a last 
flight to be optimised has been processed. If not, the 
method returns to block 308 for applying the operations 
15 described above to a next flight. Conversely, the method 
returns to block 304 waiting for a new command. 

Similar considerations apply if the program 
application performs an equivalent method, for example 

2 0 with error routines, escape functions for stopping 

execution of the program, several concurrent processes 
that execute the above described operations in parallel, 
and the like. Alternatively, a resolution of the freights 
or a unit load device used for transporting the freights 

2 5 are considered (for example by reducing the residual 
capacity of a pre-set coefficient) , no-show freights 
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and/or denied boarding (off-loaded) freights are 
estimated, an overbooking is taken into account (for 
example increasing the offered capacity by a pre-set 
coefficient) , or no by-product is calculated. 

5 

More generally, the present invention provides a 
yield management method, which is used to optimise a 
yield parameter resulting from assigning a capacity 
offered by a future instance of a service to each one of 

10 a plurality of different categories of requests competing 
for the capacity; the capacity is defined by one or more 
capacity variables. A set of historical profiles is 
stored for each one of a plurality of previous instances 
of the service; the set includes a historical profile of 

15 an historical value of each capacity variable reserved by 
each category. A probability is then assigned to each 
previous instance of the service. The method estimates a 
potential profile of a potential value of the capacity 
variable from each historical profile; the estimation is 

2 0 carried out according to a corresponding current value of 
the capacity variable reserved by the category for the 
future instance of the service and according to a 
corresponding unconstrained demand of the capacity 
variable for the category in the previous instance of the 

25 service. A historical scenario is defined for each 

previous instance of the service; the historical scenario 
includes a final potential capacity variable from each 
corresponding potential profile. The method then 
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determines an authorisation to allocate the offered 
capacity for each capacity variable of each category in 
the future instance of the service by applying a 
stochastic model to the historical scenarios according to 
5 the corresponding probabilities. 

The computational complexity of the method according 
to the invention increases linearly with the number of 
the historical scenarios and it is independent of the 

10 number of categories; this is particularly advantageous 

in a commercial application wherein the market dynamic is 
very fast, so that a relatively low number of historical 
scenarios is often enough to take planning decisions. As 
a consequence, the problem defined by the yield 

15 management method can be solved in a reasonable amount of 
time. This allows the nesting authorisations to be 
recalculated several times during the opening period, and 
particularly with a higher frequency close to the 
departure of the future instance of the flight. 

20 

The solution of the invention does not require any 
simplification, shortcut or assumption to be made. 
Particularly, all the capacity variables are taken into 
account simultaneously, the scenarios are not 
25 artificially reduced in any way, and no assumption is 

made on the (simple and joint) statistical distribution 
of the variables involved in the problem. 
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The envisaged solution derives the probabilistic 
nature of the demand directly from the reality. As a 
consequence, any kind of demand {even if particularly 
scattered or inconsistent) yields to consistent and 
5 conservative suggestions. Therefore, the method of the 
invention is accurate, robust and secure. 

The preferred embodiment of the invention described 
above offers further advantages. For example, the 

10 probability assigned to each historical scenario is 

calculated by normalising the corresponding weight, which 
is a combination of the season coefficient and the time 
distance coefficient. In this way, the previous instances 
of the flight that have happened in a similar period as 

15 their future instance and/or that are close to the 
current time are given more importance. 

Preferably, this value is calculated as a weighted 
sum of the season coefficient and the time distance 
2 0 coefficient. This approach gives more relevance to the 

dominant coefficient, and yields to a good result even if 
the number of previous instances of the flight is not 
very high. 

2 5 Similar considerations apply if a different number 

of periods is taken into considerations (such as the days 
of the week) , if the season coefficient and/or the time 
distance coefficient are calculated in a different 
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manner, if different values are used in the formula 
defining the time distance coefficient and/or the weight, 
if the time distance coefficient depends on a release of 
the previous instance of the flight (or more generally on 
5 an occurrence thereof) , and so on. Alternatively, each 
weight is assigned to a set consisting of more than one 
previous instance of the flight (such as including all 
the previous instances of the flight in a corresponding 
month) , the weight further depends on mobile events (such 
10 as Easter) , the weight is calculated as the product of 

the season coefficient and the time distance coefficient 
(giving a high weight only when both the coefficients are 
high) . 



15 The potential profiles are estimated by emphasising 

the demand when the respective category has been closed 
to sales; the emphasis values are defined in order to 
take into account both the weight and the volume in a 
symmetric manner. The suggested procedure emphasises the 

2 0 demand according to what has happened in the other 

previous instances of the flight in which the category 
was (almost) always open. 



Moreover, the result of the integration is corrected 
2 5 to ensure that no negative forecasts for the potential 
requests at the departure of the future instance of the 
flight are obtained when the demand is low. 
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Similar considerations apply if each profile 
consists of a different number of snapshots, if the 
snapshots are taken with a different frequency, if the 
potential profiles are built with a different procedure, 
5 if the opening coefficients and/or the emphasis values 
are calculated in another manner, if a different 
interpolation technique for obtaining the potential 
gradients is employed, if the potential gradients are 
integrated in a different manner, if a dummy snapshot 
10 preceding the opening of the flight is introduced for 

processing flights that have not been opened yet, and the 
like. 

Alternatively, the emphasis values are calculated 
15 taking into account only the previous instances of the 

flight in which the category was always open (assuming a 
sufficient number of corresponding historical profiles) , 
using only the historical profiles with a growing demand, 
or performing a weighted mean of the more optimistic 
20 increments of the demand in the historical profiles. 

However, the method according to the present invention 
leads itself to be carried out without any correction to 
the result of the integration, or even detecting and 
storing any denied request, so as to know the potential 
25 profiles directly (even if some adjustment is always 

suggested because it is very difficult to known all the 
denied requests for sure) . 
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The potential profiles are further reconciled to an 
admissible reference density; in this way, it is assured 
that no inconsistent forecasts are produced. 

5 The reference density is set to the closest one 

between the density of the previous snapshot in the 
potential profile and the density of the current snapshot 
in the historical profile. This value makes it possible 
to take into account what has happened in the historical 
10 profile without excessively affecting the potential 
profile . 



Similar considerations apply if another logic 
relation between the capacity variables is considered for 

15 the reconciliation of the potential profiles, if the 

reconciliation is carried out in a different manner, if a 
different reference density is employed, and so on. For 
example, in an alternative embodiment of the present 
invention the reference density is set to a mean of the 

20 above-mentioned limit values, to the density of the 

preceding snapshot in the historical profile, or to a 
pre-set value. However, the method of the invention leads 
itself to be carried out even without performing any 
reconciliation of the potential profiles. 



Each historical scenario further includes the 
potential revenue, which is estimated from the 
corresponding historical revenue; when the value is 
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unknown, the historical revenue is set to a weighted mean 
of the values known in the other previous instances of 
the flight. This is particularly advantageous when it is 
difficult, if not impossible, to detect the real revenue 
5 of each shipment in the previous instance of the flight 
at its departure; moreover, the employed procedure is 
independent of the hypothesis that the freights of each 
shipment maintain the same category from the historical 
profile to the potential profile. 

10 

Preferably, the potential unit revenue is calculated 
as a mean of the historical unit revenue and the current 
unit revenue, weighted according to the corrective 
factor. In this way, the potential unit revenue takes 
15 into account what is currently happening for the future 

instance of the flight; however, the current unit revenue 
affects the potential unit revenue only when the 
corrective factor is high (this is preferable when the 
current revenue is not particularly accurate) . 

20 

The corrective factor is a combination of the 
departure coefficient and the increment coefficient. In 
this way, the potential revenue takes into account what 
is actually happening only when this can affect the 
25 forecast. 

The corrective factor is calculated as a weighted 
sum of the departure coefficient and the increment 
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coefficient. This approach gives more relevance to the 
dominant coefficient, and yields to a good result even if 
the current revenue is not very accurate. 

5 Similar considerations apply if the potential 

revenue is estimated in a different manner, if the value 
of each accepted request is set to a mean market value 
(when its value cannot be determined in an accurate 
manner by the reservation subsystem) , if another formula 

10 is employed for combining the potential unit revenue and 
the current unit revenue, if the departure coefficient 
and/or the increment coefficient are calculated in a 
different manner, if different values are used in the 
formula defining the departure coefficient and/or the 

15 corrective factor, and so on. Alternatively, the 

corrective factor depends on other coefficients, the 
corrective factor is calculated as the product of the 
departure coefficient and the increment coefficient 
(giving a high value only when both the coefficients are 

20 high); however, the solution of the invention leads 
itself to be implemented without taking the current 
revenue into consideration, or even without estimating 
any potential revenue. 

2 5 Two nesting orders (for the weight and the volume) 

are calculated according to the mean potential revenue of 
each category. 
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Moreover, corresponding input nesting orders are 
also provided; the nesting orders are obtained updating 
the input nesting orders according to the mean potential 
revenue for the categories in which this value is 
5 available. 

Similar considerations apply if another technique is 
employed for calculating the nesting orders, or if the 
input nesting order are updated in a different manner. 
10 Alternatively, no input nesting orders are provided, or 
the nesting orders are not calculated by the system (but 
they are always set by the space controller) . 

The nesting portion and the out of nesting portion 
15 of the target function, which are combined using the 
aggressiveness parameter, allow the result of the 
optimisation to be tuned according to different attitudes 
to the risk; particularly, a more aggressive behaviour 
gives a higher revenue in the worst case (wherein the 
2 0 requests for the most valuable freights are the last 

ones) , and a lower revenue in the best case (wherein the 
requests for the most valuable freights are the first 
ones) . 

25 Advantageously, the out of nesting portion of the 

target function is calculated with two independent 
components for the weight and for the volume. In this 
way, the constraints for the worst case are slightly 
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relaxed with respect to a real out of nesting policy, 
wherein the out of nesting authorisations must be met by 
both the weight and the volume. Conversely, in the 
proposed formula the two capacity variables are 
5 considered one at the time, so that each request may be 
accepted even if it meets the out of nesting 
authorisations only for the weight and not for the volume 
(or vice-versa) . 

10 Similar considerations apply if the stochastic model 

is defined in a different manner, if another target 
function is envisaged, if the different components of the 
out of nesting portion and/or the two portions of the 
target function are combined in another way, and the 

15 like. However, the solution of the invention leads itself 
to be carried out even with the out of nesting portion 
that is calculated simultaneously for the weight and the 
volume, or with the target function that is defined with 
a single policy (either in nesting or out of nesting) . 

20 

Preferably, the space controller inputs a user 
scenario with the respective probability; in this way, 
any personal feeling of the future may be introduced into 
the stochastic model in a very simple manner. However, 
2 5 the solution according to the present invention leads 

itself to be implemented even without the possibility of 
imputing any further scenario. 



GB920010034US1 



51 



The proposed solution is particularly advantageous 
in a cargo flight, wherein the capacity is defined by two 
independent variables (weight and volume) , and the 
revenue of the old instances of the flight is not 
5 immediately available. Alternatively, the solution of the 
invention is used to optimise another yield parameter 
resulting from assigning an offered capacity defined by 
one or more different variables; moreover, other 
applications of the yield management method are 
10 contemplated, for example in a passenger flight, in a 

different transportation system (such as a train) , in a 
system for controlling booking of rooms in a hotel, and 
the like. 

15 Advantageously, the solution according to the 

present invention is implemented with a program 
application, which consists of a computer program 
(software) provided on CD-ROM. 

2 0 Alternatively, the program application is provided 

on floppy-disk, is pre-loaded onto the hard-disk, or is 
stored on any other computer readable medium, is sent to 
the authorisation subsystem through the INTERNET, is 
broadcast, or more generally is provided in any other 

2 5 form directly loadable into the working memory of the 

authorisation subsystem. However, the method according to 
the present invention leads itself to be carried out even 
with a hardware structure, for example integrated in a 



GB920010034US1 



52 



chip of semiconductor material . 

Naturally, in order to satisfy local and specific 
requirements, a person skilled in the art may apply to 
the solution described above many modifications and 
alterations all of which, however, are included within 
the scope of protection of the invention as defined by 
the following claims . 



